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ABSTRACT 

For all available pseudo-Schwarzschild potentials, we provide a non-self-similar model 
of coupled accretion-outflow system in connection to the Quasi Periodic Oscillation 
(QPO) of the black hole powered galactic microquasars and the emergence of barionic 
jets out of these objects. We use the vertically integrated 1.5 dimensional model to 
describe the disc structure where the equations of motion are written on the equatorial 
plane of the central accretor, assuming the flow to be in hydrostatic equilibrium in the 
transverse direction. First we formulate and solve the equations governing axisymmet- 
rically rotating, advective, multi-transonic black hole accretion which may contain the 
Rankine Hugoniot Shock Waves (RHSW), and then we calculate the associated QPO 
frequencies vqpo in terms of relevant accretion parameters. We then argue that the 
post-shock region for such flows may serve as an efficient source of outflow generation; 
and we calculate, for same set of accretion parameters used to calculate the vqpos, 
what fraction of the accreting material, denoted by R m , is being blown as shock gen- 
erated outflow. In this way we theoretically study the relation between vqpos and 
R m s for galactic microquasars and compare our theoretically obtained result with 
observational data. 
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1 INTRODUCTION 

It has been established in recent years that in order to satisfy the inner boundary conditions imposed by the event horizon, 
accretion onto black holes should exhibit transonic properties in general; which further indicates that formation of shock 
waves are possible in astrophysical fluid flows onto galactic and extra-galactic black holes. One also expects that shock for- 
mation in black hole accretion might be a general phenomena because shock waves in rotating and non-rotating flows are 
convincingly able to provide an important and efficient mechanism for conversion of significant amount of the gravitational 
energy (available from deep potential wells created by these massive compact accretors) into radiation by randomizing the 
directed infall motion of the accreting fluid. Hence shocks possibly play an important role in governing the overall dynamical 
and radiative processes taking place in astrophysical fluid and plasma accreting onto black holes. Thus the study of steady, 
standing, stationary shock waves produced in black hole accretion has acquired a very important status in recent years and 
it is now believed that shocks may be an important ingredient in an accreting black hole system in general (for an extended 
list of literatures on shock formation in black hole accretion, see Das 2002, D02 hereafter, and references therein). Hot, dense 
and exo-entropic post-shock regions in advective accretion disks are also used as a powerful tool in understanding the spec- 
tral properties of black hole candidates (Shrader & Titarchuk 1998, and references therein) and in theoretically explaining a 
number of diverse phenomena, including the generation mechanism for high frequency Quasi Periodic Oscillations (QPOs) in 
general (Titarchuk, Lapidus & Muslimov 1998, and references therein). A number of observational evidences are also present 
which are in close agreement with the theoretical predictions obtained from shocked accretion model (Rutledge et al. 1999; 
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Muno, Morgan & Remillard 1999; Webb & Malkan 2000; Rao, Yadav & Paul 2000; Smith, Heindl & Swank 2001). In their 
attempt to explain some of the observational features of galactic microquasars, Chakrabarti & Manickam (Chakrabarti & 
Manickam 2000) proposed that the intermediate frequency (z/~l — WHZ) and high frequency QPO of the black hole candidate 
(harbored by GRS 1915+105) may occur due to oscillations of shocks in accretion disc. In a recent work, Das (2003, hereafter 
D03), showed that the QPOs in galactic sources are, indeed, regulated by shocked accretion flow, and analytically calculated 
the QPO frequency as a function of fundamental accretion parameters. It is now a well established fact that quasars and 
microquasars suffer mass loss through outflows and jets (Mirabel & Rodriguez 1999; Ferrari 1998; Begelman, Blandford & Rees 
1984). These galactic and extra-galactic jet sources are commonly believed to harbor accreting compact objects at their hearts 
as the prime movers for almost all non stellar energetic activities around them including the production of bipolar outflows 
and relativistic jets. Unlike normal stellar bodies, compact objects do not have their own physical atmosphere from where 
matter could be ripped off as winds, hence outflows from the vicinity of these prime movers have to be generated only from 
the accreting material. So instead of separately investigating the jets and accretion processes as two disjoint issues around the 
dynamical center of the galactic and extra-galactic jet sources, it is very necessary to study these two phenomena within the 
same framework and any consistent theoretical model for jet production should explore the outflow formation only from the 
knowledge of accretion parameters. Also to be noted that while self-similar models are a valuable first step, they can never be 
the full answer, and indeed any model which works equally well at all radii is fairly unsatisfactory to prove its viability. Thus 
the preferred model for jet formation must be one which is able to select the specific region of jet formation. Motivated by 
the above mentioned arguments, Das (Das 1998, D98 hereafter, Das 2001) and Das & Chakrabarti (Das & Chakrabarti 1999, 
DC hereafter) proposed a non self-similar analytical model capable of self-consistently exploring the hydrodynamic origin of 
accretion powered jets/outflows emanating out from galactic and extra-galactic sources. D98 and DC computed mass outflow 
rates from accretion disks around compact objects, such as neutron stars and black holes. These computations were done 
using combinations of exact transonic inflow and outflow solutions which may form standing shock waves. Assuming that 
the bulk of the outflow is from the shock generated effective boundary layers of these objects, they found that the ratio of 
the outflow rate to the inflow rate varies anywhere from a few percent to even close to a hundred percent (i.e., close to disk 
evacuation case) depending on the initial parameters of the disk, the degree of compression of matter near the centrifugal 
barrier, and the polytropic index of the flow. Also the exact location (distance measured from the central accretor) of the jet 
launching zone has been successfully pointed out in their work. 

Since D98, D03 and DC share a common underlying phenomena, the formation of steady, standing Rankine Hugoniot shock 
waves (RHSW hereafter) in a thin, rotating, axisymmetric, inviscid steady accretion flow around black holes, we believe that it 
is worth investigating whether theoretical calculations allows one to correlate the frequencies of QPO of a black hole candidate 
with the amount of barionic content of jets emanating from the vicinity of that black hole sitting at the heart of any galactic 
microquasar, and if so, whether such theoretical correlation may further be supported by present day observational evidences. 
Our approach in this paper is precisely this. First we formulate and solve the equations governing axisymmetrically rotating, 
multi-transonic inviscid black hole accretion which may contain RHSW and then we calculate the associated QPO frequencies 
vqpo in terms of relevant accretion parameters. We then argue that the post-shock region for such flows may serve as an 
efficient source of outflow generation and we calculate, for the same set of accretion parameters used to calculate the vqpos, 
what fraction of the accreting material, denoted by Rm, is being blown as shock generated wind. Ultimately we provide the 
correlation among vqpos and RmS and check the feasibility of our theoretically obtained result against observational data. 



2 INFLOW EQUATIONS, MULTIPLICITY OF SONIC POINTS AND SHOCK FORMATION 

Rigorous investigation of the complete general relativistic multi-transonic black hole accretion disc and wind is believed to 
be extremely complicated. At the same time it is understood that as relativistic effects play an important role in the regions 
close to the accreting black hole (where most of the gravitational potential energy is released), purely Newtonian gravitational 
potential (in the form $jv = — GM ^ H , where Mbh is the mass of the compact accretor) cannot be a realistic choice to 
describe transonic black hole accretion in general. To compromise between the ease of handling of a Newtonian description of 
gravity and the realistic situations described by complicated general relativistic calculations, a series of 'modified' Newtonian 
potentials have been introduced to describe the general relativistic effects that are most important for accretion disk structure 
around Schwarzschild and Kerr black holes. Introduction of such potentials allows one to investigate the complicated physical 
processes taking place in disc accretion in a semi-Newtonian framework by avoiding pure general relativistic calculations so 
that most of the features of spacetime around a compact object are retained and some crucial properties of the analogous 
relativistic solutions of disc structure could be reproduced with high accuracy. Hence, those potentials might be designated 
as 'pseudo-Kerr' or 'pseudo- Schwarzschild' potentials, depending on whether they are used to mimic the space time around 
a rapidly rotating or non rotating/ slowly rotating (Kerr parameter a ~ 0) black hole respectively. 

It is important to note that as long as one is not interested in astrophysical processes extremely close (within 1 — 2 r g ) 
to a black hole horizon, one may safely use these black hole potentials to study accretion on to a Schwarzschild black hole 
with the advantage that use of these potentials would simplify calculations by allowing one to use some basic features of flat 
geometry (additivity of energy or de-coupling of various energy components etc.) which is not possible for calculations in a 
purely Schwarzschild metric. Also, one can study more complex many body problems such as accretion from an ensemble of 
companions or overall efficiency of accretion onto an ensemble of black holes in a galaxy or for studying numerical hydrodynamic 
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accretion flows around a black hole etc. as simply as can be done in a Newtonian framework, but with far better accuracy. 
However, one should be careful in using these potentials to study shocked black hole accretion because of the fact that none 
of the potentials discussed here are 'exact' in a sense that they are not directly derivable from the Einstein equations. These 
potentials could only be used to obtain more accurate correction terms over and above the pure Newtonian results and any 
'radically' new results obtained using these potentials should be cross-checked very carefully with the exact general relativistic 
theory. 

Our calculations in this paper are based on four such pseudo Schwarzschild potentials. Explicit forms of those four potentials 
are the following (see Das & Sarkar 2001 and D02 references therein for detail discussion about various properties of these 
potentials) : 



*i(r) = - 



2(r-l) 



*a(r) = 



j_ 

' 2r 



_3_ 

2r 



$ 3 (r) = - 




where r is the radial co-ordinate scaled in units of Schwarzschild radius. Hereafter, we will define the Schwarzschild radius r g 

as 

2GMbh 

(where Mbh is the mass of the black hole, G is universal gravitational constant and c is velocity of light in vacuum) so that 
the marginally bound circular orbit r& and the last stable circular orbit r s take the values 2r 3 and 3r s respectively for a typical 
Schwarzschild black hole. Also, total mechanical energy per unit mass on r s (sometimes called 'efficiency' e) may be computed 
as —0.057 for this case. We will use a simplified geometric unit throughout this paper where radial distance r is scaled in 
units of r g , radial dynamical velocity u and polytropic sound speed a of the flow is scaled in units of c (the velocity of light 
in vacuum), mass m is scaled in units of Mbh and all other derived quantities would be scaled accordingly. For simplicity, 
we will use G = c = 1. 

Among the above potentials, 3>i(r) was introduced by Paczyhski and Wiita (1980) which accurately reproduces the positions 
of r 3 and rv Also the Keplerian distribution of angular momentum obtained using this potential is exactly same as that 
obtained in pure Schwarzschild geometry. 3>2(r) was proposed by Nowak and Wagoner (1991) to approximate some of the 
dominant relativistic effects of the accreting black hole (slowly rotating or non-rotating) via a modified Newtonian potential. 
It has the correct form of r s as well as it produces the best approximation for the value of the angular velocity fi s (as measured 
at infinity) at r 3 and the radial epicyclic frequency k (for r > r s ). <£>3(r) and $4(r) were proposed by Artemova, Bjornsson & 
Novikov 1996, ABN hereafter, to produce exactly the same value of the free-fall acceleration of a test particle at a given value 
of r as is obtained for a test particle at rest with respect to the Schwarzschild reference frame ($3) and to produce the value 
of the free fall acceleration that is equal to the value of the covariant component of the three dimensional free-fall acceleration 
vector of a test particle that is at rest in the Schwarzschild reference frame ($4) respectively. Hereafter, we will denote any 
ith potential as $i(r) where {i = 1,2,3,4} corresponds to {$1 (r), $2(7"), $3(»"), $4(7*)} respectively. 

Following standard literature, we consider a thin, rotating, axisymmetric, inviscid steady flow in hydrostatic equilibrium in 
transverse direction. The assumption of hydrostatic equilibrium is justified for a thin flow because for such a flow, the infall 
time scale is expected to exceed the local sound crossing time scale in the direction transverse to the flow. The flow is also 
assumed to possess considerably large radial velocity which makes the flow 'advective' (see D02, and references therein). The 
complete solutions of such a system (which may allow RHSW to form) require the dimensionless equations for conserved 
specific energy £ and angular momentum A of the accreting material and the mass conservation equations supplied by the 
transonic conditions at the sonic points and the Rankine Hugoniot conditions at the shock. The local half-thickness, hi (r) of 
the disc for any $i can be obtained by balancing the gravitational force by pressure gradient and can be expressed as: 

fti(r) = oVr/(7^) (2), 

<3^ is the derivative of any ith potential with respect to the radial co-ordinate. For a non-viscous flow obeying the polytropic 
equation of state p — Kp~' (K is a measure of the specific entropy of the flow), integration of radial momentum equation: 

du 1 dp d r^eff, ,n 

u— + --f + — m (r)\ = 
dr pdr dr 1 J 

where $^ (r) is the iih 'effective' potential which is the summation of the gravitational potential and the centrifugal potential 
for matter accreting under the influence of ith pseudo potential: 

leads to the following energy conservation equation in steady state: 

£ = ^ + ^T + ^+^W=0; (3a) 
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Similarly, the continuity equation: 

\uprhi(r)] = 

dr 

can be integrated to obtain the barion number conservation equation as: 

M = ^u e a e p e r'i . (3b) 

■ (-*-) (-*-) 

Following standard literature, one can define the entropy accretion rate M = MK V t- 1 / 7V / which undergoes a discon- 
tinuous transition at the shock location r s h where local turbulence generates entropy to increase M for post-shock flows. For 
our purpose, explicit expression for M can be obtained as: 

M = y^iteai^r-i . (3c) 

In Eqs. (3a-3c), the subscript e indicates the values measured on the equatorial plane of the disk; however, we will drop e 
hereafter if no confusion arises in doing so. 

At this point, we would like to discuss in detail about some of the assumptions made in this work to study the disc structure 
and dynamics. Firstly, our calculations are essentially based on inviscid accretion whereas in reality one may expect that 
viscosity is present in accretion disc around black holes. The exact expression of A in eq. (8a) may, thus be A(r) (which means 
that the angular momentum would be a function of the radial distance) and along with the radial momentum equation and 
continuity equation, two other equations, namely, equations for azimuthal angular momentum and the equation for heating 
and cooling factors describing the radiative properties of the flow, should be used to bring the whole picture into focus. 
However, we try to justify our assumption of inviscid flow in the following way: 

We will mainly concentrate in a region of accretion disc close to the central accretor (a few Schwarzschild radii away from 
the event horizon). This is because, except for a few unusual situations, for all reasonable initial boundary conditions for 
flow in every 3> 4 , shocks are likely to form in this length scale (see D02). Close to the black hole at that length scale, the 
radial velocity of matter becomes enormously large due to extremely strong gravitational attraction of the black hole, hence 
infall time scale becomes much smaller compared to the viscous time scale for all practical purposes. So, one can treat the 
angular momentum to be practically constant irrespective of the nature of viscosity in accretion disc and our assumption of 
inviscid flow may not be quite unjustified. Nevertheless, one limitation of our model is to treat the flow to be inviscid even 
at far field position, i.e, at a long distance away from the accretor where the viscous dissipation of angular momentum may 
be taken into account. We thus, were unable to prescribe any analytical means using which we could smoothly join the far 
field high viscous flow with near field weakly viscous / practically inviscid flow. However, even thirty years after the discovery 
of standard accretion disc theory (Shakura & Sunyaev 1973), exact modeling of viscous multi-transonic black hole accretion, 
including proper heating and cooling mechanism is quite an arduous task, and we do not attempt it in this paper. However, 
our qualitative calculations show that the introduction of viscosity via a radius dependent power law distribution for angular 
momentum only weakness the strength of the centrifugal barrier and pushes the shock location closer to the event horizon, 
keeping the overall basic physics concerning the shock dynamics and related issues unaltered; details of this work is in progress 
and will be discussed elsewhere. 

One can simultaneously solve Eqs. (3a-3c) for any particular <3>i and for a particular set of values of {£, A, 7}. Hereafter we 
will use the notation [Pi] for a set of values of {£,A,7} for any particular $ 4 . 

For a particular [Pi], it is now quite straight-forward to derive the space gradient of dynamical flow velocity for flow in 

any particular ith black hole potential <!>;(r) as: 



(^ + *i(r))-^l7t - 



= ^ 



2a 2 
«(7+l) 

Since the flow is assumed to be smooth everywhere, if the denominator of eqn. (4a) vanishes at any radial distance r, the 
numerator must also vanish there to maintain the continuity of the flow. One therefore arrives at the so called 'sonic point 
(alternately, the 'critical point') conditions' by simultaneously making the numerator and denominator of eqn. (4a) equal to 
zero. The sonic point conditions then can be expressed as: 



1 + 7 i 
—z — u s 



3$;(r) + r$"(r) 



(46) 



where the subscript s indicates that the quantities are to be measured at the sonic point (s) and $/' represents the derivative 
of For a fixed [Pi] and $ 4 , one can solve the following polynomial of r to obtain the sonic point(s) of the flow where the 
subscript s refers to the quantities measured at the sonic points: 



2r 2 y 7 2 



*iW +l<$> t (r) ( A 2 +r 3 $ I (r) 



0. (4c) 
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Fig. 1: Classification of parameter space for multi-transonic accretion and wind. Shock formation region is shown for accretion only and 
not for wind. See text for details. 

Similarly, the value of (^) ; at its corresponding sonic point (s) r\ can be obtained by solving the following equation: 



47 / du 



7+1 \ dr 

Cm 27 



2u s 



7+1 



3 | Mr) 



(du\ 
\dr J s a 



6(7-1) /$i(r) 



6 (2 7 - 1) 

(l + 7) 2 V*iWy 7(7+l)U^)y 7 2 (7+l) 2 



+$, 



3A" 

^4 



(4d) 



Where the subscript "s, i" indicates that the corresponding quantities for any iih potential is being measured at its corre- 
sponding sonic point (s) and 3^ (r) — —§pt-^-- 

For all <E>i's, we find a significant region of parameter space spanned by [Pi] which allows the multiplicity of sonic points for 
accretion as well as for wind where two real physical inner and outer (with respect to the black hole location) X type sonic 
points ri n and r out encompass one O type unphysical middle sonic point r mi£ j in between them. If shock forms in accretion 
(in this work we will not study the shock formation in wind), then [Pi]s responsible for shock formation must be somewhere 
from the region for which three sonic points will form in accretion only, though not all [Pi]s in the region of multi-transonic 
accretion will allow shock transition (see subsequent discussions and Fig. 1). Using Eqs. (2,3a-3c), one can combine the three 
standard RH conditions (Landau & Lifshitz 1959) for vertically integrated pressure and density to derive the following relation 
which is valid only at the shock location: 



(1-7) 



P-M- 
M 



log L r 



'(ki+th) 



9(1 + 9 - RcomvT 1 + (1 + 6)" 



o, 



(5) 



where £(ki+th) is the total specific thermal plus mechanical energy of the accreting fluid: £(ki+th) = £ — {^pi + 

Rcomp and P are the density compression and entropy enhancement ratio respectively, defined as Rcomp = (p+/p-) an d 
P = (.M + /.M-) respectively; 6 = 1 — r' 1-7 - 1 and Y = (3R comp , "+" and "_" refer to the post- and pre-shock quantities. The 
shock strength Si (ratio of the pre- to post-shock Mach number of the flow) can be calculated as: 



Si — Rcomp (1 + 0) • 



(6) 



One can simultaneously solve Eqs. (3a-3c,5-6) to find out the shock location r s h along with any sonic or shock quantity as a 
function of [Pi] and can identify the regions of parameter space (spanned by £, A and 7) responsible for shock formation for 
any 5>i(r). However, one can show that as far as the shock formation in multi-transonic accretion flows around black holes 
are concerned, the Paczyhski and Wiita potential is expected to be the best approximation of complete general relativistic 
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solutions (D02), and hence, although we will provide a general solution scheme for shock induced outflow for all pseudo 
potentials, we will present our main results for accretion in $i(r) only. 

In Fig. 1, we classify the parameter space for multi-transonic accretion and wind for flows in 3>i(r). The specific energy £ is 
plotted along the Y axis and the specific angular momentum A is plotted along the X axis. In the region bounded by PQR, 
three sonic points are formed in accretion, while three sonic points in wind are formed in the region bounded by PRS S. 
The shaded wedge PTR represents the region of parameter space for which multi-transonic accretion flow will have a steady, 
standing RHSW. Any region which falls outside PQS S, will produce a mono-transonic accretion, so to say. Although the 
figure is drawn for ultra-relativistic * flow, one can easily explore multi-transonic shocked accretion flow with higher values 
of 7 (4/3 < 7 < 5/3). For any [Pi] £ PTR, one can calculate the shock location r sh and the shock compression ratio R C om P - 
Following D03, one can calculate the QPO frequency vqpo as: 

vqpo = (7) 

Rcompr sh 

The dimensionless constant A in the above equation can be determined by scaling the calculated QPO frequencies with the 
highest possible observed vqpo for a particular astrophysical source. 



3 GENERATION OF SHOCK-INDUCED OUTFLOW AND ITS GOVERNING EQUATIONS 

The simplicity of black holes lie in the fact that they do not have atmospheres. But the advective accretion discs surrounding 
them have, and similar method as employed in stellar atmospheres should be applicable to the disks. Our approach in this 
section is precisely this. We first explain how the post shock region is expected to produce outflows; in doing so we will follow 
more or less the same arguments used by D98 and DC. Then we will formulate the equations governing such outflows in a 
general format so that we will incorporate all pseudo potentials. Finally we will simultaneously solve the equations governing 
the inflow and outflow to find out what fraction of the accreting material, denoted by Rm, is being blown as wind and for 
the same set of accretion parameters, we will correlate that fractional amount with the frequencies of the QPO. Due to the 
fact that close to the black hole the radial component of the infall velocity of accreting material would be enormously high, 
viscous time scale would be much longer than the infall time scale and a rotating inflow entering into a black hole will have 
almost constant specific angular momentum close to the black hole for any moderate viscous stress. This almost constant 
angular momentum produces a very strong centrifugal force which increases much faster compared to the gravitational force 
and becomes comparable at some specific radial distance. Here, (actually, a little farther out, due to thermal pressure) matter 
starts piling up and produces the centrifugal pressure supported boundary layer (CENBOL). Further close to the black 
hole, the gravity always wins and matter enters the horizon supersonically after passing through a sonic point. Formation of 
CENBOL may be attributed to the shock formation in accreting fluid. In CENBOL region the post-shock flow becomes hotter 
and denser and for all practical purposes behaves as the stellar atmosphere so far as the formation of outflows are concerned. 
A part of the hot and dense shock-compressed inflowing material is then 'squirt' as outflow from the CENBOL. Subsonic 
outflows originating from CENBOL would pass through outflow sonic points and reach far distances as in wind solution. 
It is to be noted here that the generation of outflow from CENBOL in this work is a rational assumption. The exact analytical 
calculation describing the change of linear momentum of the accreting material in a direction perpendicular to the plane of the 
disc is beyond the frame work of the 1.5 dimensional disc model used in this work where the explicit variation of dynamical 
variables along the Z axis is not amenable to analytical treatment. It can be stated that the enormous post shock thermal 
pressure in CENBOL is, in reality, capable of providing a substantial amount of 'hard push' to the accreting material against 
the gravitational attraction of the black hole and this 'thermal kick' may play an important role in re-distributing the linear 
momentum of the inflow and generates a non zero component along the Z direction. In other words, thermal pressure at 
CENBOL, being anisotropic in nature, may deflect a part of inflow along the perpendicular to the equatorial plane of the 
disc, exact mechanism of which we could not formulate analytically. Thus in our model, the CENBOL successfuly connects 
the two stationary quasi-one dimensional solutions, namely, the 1.5 dimensions vertically averaged inflow and the outflow, to 
produce a coupled disc-outflow system. 

However, we do believe that the CENBOL region is generic (unless viscosity is exorbitantly high which removes angular 
momentum of the flow almost completely), because all the self-consistent solutions of the governing equations show that 
angular momentum is almost constant close to the black hole even though the matter start with a large angular momentum 
at a large (few million) Schwarzschild radius away. As we explained, the viscous time scale is always large compared with the 
infall time scale, and thus angular momentum remains almost constant and produces the CENBOL. We are not assuming that 
usual Keplerian disks (farther away) do not form. We think that while forming this Keplerian disk, the angular momentum 
has been re-distributed and excess has been absorbed by the companion as in the usual scenario. What is new here is the 
realization that the angular momentum distribution must deviate from a Keplerian distribution close to the black hole and 
therefore produces a CENBOL closer to the black hole. Such a picture of the inflow with Keplerian flow in the high viscosity 
region of equatorial plane and low viscosity flow away from it is not due to convenience or necessity, but it arises out of 

* By the term 'ultra-relativistic' and 'purely non-relativistic' we mean a flow with 7 = § and 7 = f respectively, according to the 
terminology used in Frank et. al (1992). 
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Fig. 2: Multi-component combined flow geometry in 3-Dimension 

the actual solutions of the viscous transonic flow equations (see Chakrabarti & Titarchuk 1995 and references therein). This 
has not only changed our view of the state changes of the black hole spectra (hard/soft), this has also given us clues to the 
power-law tail of the soft state etc. through the bulk motion Comptonization etc. The interesting feature is that the same 
CENBOL may be responsible for the outflow as well. 

There are two surfaces of utmost importance in flows with angular momentum. One is the 'funnel wall' where the effective 
potential $e//( r ) vanishes. In the case of a purely rotating flow, this is the 'zero pressure' surface. Flows can not enter inside 
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the funnel wall because the pressure would be negative. The other surface is called the 'centrifugal barrier'. This is the surface 
where the radial pressure gradient of a purely rotating flow vanishes and is located outside the funnel wall simply because 
the flow pressure is higher than zero on this surface. Flow with inertial pressure easily crosses this 'barrier' and either enters 
into a black hole or flows out as winds depending on its initial parameters. It is observed that (DC) the outflow generally 
hugs the 'funnel wall' and goes out in between these two surfaces. Here it is important to note that our assumption of thin 
inflow is for the sake of computation of the thermodynamic quantities only, but the flow itself need not be physically thin. 
Also the funnel wall and the centrifugal barrier are purely geometric surfaces, and they exist anyway and the outflow could 
be supported even by ambient medium which may not necessarily be a part of the disk itself. It is true that the surfaces could 
have been constructed more accurately if the viscous and other dissipative forces could also be used to obtain this. But the 
error cannot be high, since the radial motion in the outflow is very small compared to the azimuthal motion at least up to 
the sonic point. 

Qualitatively, one might attempt to 'visualize' how the combined accretion-outflow system along with the central accretor 
would 'look like' in reality. In figure 2., we attempt to illustrate the 3-D geometry of coupled disk-outflow system according 
to our model. B is the accreting Schwarzschild black hole while C's represent the hot and dense CENBOL region. D(K) 
and D(SK) represent the thin Keplerian part and puffed- up sub-Keplerian part of the advective accretion disk respectively 
(diagram not in scale). Due to the axisymmetry assumption in accretion, two oppositely directed jet are expelled from the 
close vicinity of B. The inner and the outer surfaces of the outflow are the funnel wall and centrifugal barrier respectively as 
explained above, see eq. (9a-9b) in next paragraph. 

In ordinary stellar mass loss computations, the outflow is assumed to be isothermal till the sonic point (Tarafdar 1988, and 
references therein). This assumption is probably justified, since copious photons from the stellar atmosphere deposit momenta 
on the slowly outgoing and expanding outflow and possibly make the flow close to isothermal. This need not be the case for 
outflows from compact sources. Centrifugal pressure supported boundary layers close to the black hole are very hot (close to 
the virial temperature) and most of the photons emitted may be swallowed by the black holes themselves instead of coming 
out of the region and depositing momentum onto the outflow. Thus, the outflows could be cooler than isothermal flows. In 
our work, we choose polytropic outflows with the same energy as the inflow (i.e., no energy dissipation between the inflow and 
outflow) but with a different polytropic index 7 < 7. We have approximated 70 as a free parameter though in reality j is 
directly related to the heating and cooling processes of the outflow. However, investigation of such detailed radiative processes 
are beyond the scope of this work. When Mi n is high, heating of outflow by photon momentum deposition is higher, and 
therefore y —> 1 letting the flow to approach towards its isothermal limit. We also assume that very little viscosity is present 
in the flow except at the place where the shock forms, so that the specific angular momentum A is constant in both inflows 
and outflows close to the black hole. That viscous time scales are longer compared to the inflow time scale, may be a good 
assumption in the disk, but it may not be a very good assumption for the outflows which are slow prior to the acceleration and 
are therefore, prone to viscous transport of angular momentum. Such detailed study has not been attempted here particularly 
because we know very little about the viscous processes taking place in the pre-jet flow. Therefore, we concentrate only those 
cases where the specific angular momentum is roughly constant when inflowing matter becomes part of the outflow. At the 
shock, entropy is generated and hence the outflow is of higher entropy for the same specific energy. 
For outflow, the energy conservation, mass conservation and the entropy conservation equations can be written as: 

„,2 ~2 \2 



2l 

7c -1 ' 2ri 



£ = - + - T + (8a) 



M = pvA(r) (86) 
M = ay^vA(r) (8c) 

where v is the outflow velocity. 

The difference between eq. (3a-3c) and (8a-8c) is that, presently, the rotational energy term contains 

Mr) = m+m, 

as the mean axial distance of the flow. The expression of TZ(r), the local radius of the centrifugal barrier comes from balancing 
the centrifugal force with the gravity: 

«<"-(&)* <9 "' 

And the expression for R(r), the local radius of the funnel wall, comes from vanishing of total effective potential, i.e.: 



*(r) = M*P (9b) 



so we can write the mean axial distance r m as: 



(9c) 
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where 4>i(r) — — &i(r). Here, A(r) is the area between the centrifugal barrier and the funnel wall. This is computed with 
the assumption that the outflow is external pressure supported, i.e., the centrifugal barrier is in pressure balanced with the 
ambient medium. Matter, if pushed hard enough, can cross centrifugal barrier in black hole accretion. An outward thermal 
force (such as provided by the enormous post shock temperature) in between the funnel wall and the centrifugal barrier causes 
the flow to come out. Thus the cross section of the outflow is, 



1 \l *;(r) 

Like accretion, the velocity gradient of the outflow can be computed as: 



(-) = 

\drJi 



J^H i 



2 



V '"« (r) 



for which the sonic point conditions comes out to be: 
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One can solve the following equation to obtain the outflow sonic point (s) 
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and the outflow velocity gradient at the sonic point(s) can be obtained by solving the following equation: 
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The mass outflow rate R m can now be defined as: 

^ = ^ = *(£,A, 7 ,7o) (11) 

where "J/ has some complicated non-linear non-linear functional form which can not be evaluated analytically. 
The following procedure is adopted to obtain a complete solution of the coupled accretion-outflow system. For any [Pi] £ PTR 
(see Fig. 1), suppose that matter first enters through the outer sonic point and passes through a shock. At the shock, part of 
the incoming matter, having higher entropy density is likely to return back as winds through a sonic point, other than the one 
it just entered. That is because, since the outflow would be heated by photons, and thus have a smaller polytropic constant, 
the flow would leave the system through an outer sonic point different from that of the incoming solution. Thus a combination 
of topologies, one from the region PTR, and the other corresponding to the flow with same (£j,Aj) but now with 7 = j a , 
are required to obtain a full solution. Thus finding a complete self-consistent solution boils down to finding the outer sonic 
point of the outflow and the mass flux through it. A supply of parameters £ , A, 7 and 7 make a self-consistent computation 
of Rm possible. We obtain the inflow sonic point by solving Eq. (4c). From Eq. (41-4d), using the fourth order Runge Kutta 



method, u(r), a(r) and the inflow Mach number 



-'('•) 



are computed along with the inflow from the inflow sonic point till the 



a(r) 

position where the shock forms. The shock location is calculated by solving Eq. (5). Various shock parameters (i.e., density, 
pressure etc at the shock surface) are then computed self-consistently. For outflow, with the known value of £ ,A and 7 , it is 
easy to compute the location of the outflow sonic point from Eq. (10c). At the outflow sonic point, the outflow velocity v\ 
and polytropic sound velocity a\ is computed from Eq. (10b). Using Eq. (10a) and (lOd), (^r) and (^) is computed as was 
done for the inflow. Runge -Kutta method is then employed to integrate from the outflow sonic point towards the black hole 
to find out the outflow velocity and density at the shock location for a given accretion rate and M ut is calculated from Eq. 
(8b). The mass outflow rate R M is then computed using Eq.(ll) and the corresponding vqpo is computed using Eq. (7). 
As $1 is the best available potential to mimic the general relativistic results for multi-transonic shocked accretion, we will 
calculate vqpo and R M only for flows in $1 by setting: 



®i(r),®i(r),$>i (r),$i (r), 



$i(r),<E' 1 (r),$ 1 (r),$! (r), 



while simultaneously solving equations (2-11), using the procedure described above. One may ask the question that if we 
plan to explore the vqpo — Rm correlation only for <J?i(r), why then we formulate all the equations for all <E>i(r)s in 
general instead of < 3>i(r), which would be a relatively less complicated procedure. Our answer to this question is mainly 
two-fold. Firstly, we wanted to present a generalized formulation for all $i(r)s to get convinced that such correlation can 
be theoretically investigated for any kind of black hole potential and to assure that our model is not just an artifact of 
a particular type of potential only. Secondly, of course there are possibilities that in future someone may come up with 
a pseudo-Schwarzschild potential better than $i(r), which will be the best approximation for complete general relativistic 
investigation of multi-transonic shocked flow. In such case, if one can in advance formulate a generalized model for coupled 
accretion- wind system for any arbitrary 3>i(r), exactly what we have done in this paper, then that generalized model will 
be able to readily accommodate that new $(r) without having any significant change in the fundamental structure of the 
formulation and solution scheme of the model and we need not have to worry about providing any new scheme exclusively 
valid only for that new potential, if any. Thus we believe that our general inflow-outflow model presented in this paper, which 
bears the same philosophy as provided by D98 and DC but treats the problem in a much more general way, will be quite 
useful for any further improvement in this field and our theoretical correlation study can also be extended easily for any 
such new and more accurate potential that might be discovered in future. It is to be noted that we cannot accelerate any 
matter more than the sound speed (or, the rotational velocity) at the initial injection point so we do not claim to explain 
superluminal motions. We have done what is hydrodynamically possible. If we follow our path for MHD flows, we believe that 
we could reach superluminal motions also. Also, we could not do a thorough job on angular momentum distribution inside 
the jet. What we understand here is that in case, some form of viscosity is working (say radiative viscosity) inside the jet to 
make the distribution of the angular momentum as a power-law, then the resulting density will go up towards the edge of the 
jet (a 'hollow' jet as in Hawley, Smarr & Wilson 1984, for example). Thus, since maximum amount of matter is taking away 
maximum amount of angular momentum (unlike in the disc), the average angular momentum need not be different from that 
of the CENBOL region. But in case, the average angular momentum were lower, then we find that the outflow rate was also 
lower. 

In Fig. 3, we represent our theoretical results regarding the correlation of vqpo with Rm- We define the normalized QPO 
frequency fnorm in the following way: 

fnorm = (12) 



V, 



max 



where v — vqpo is the theoretically calculated QPO frequency obtained from eqn. (7), with properly scaled C and v max 
is the highest possible observed frequency for any particular astrophysical source. This normalization is necessary because 
exact value of C in the equation cannot be estimated analytically. We plot fnorm along the X axis and the mass outflow rate 
scaled in solar mass per year is plotted along the Y axis. We calculate the vqpo as well R m for such shocked ultra-relativistic 
accretion on to a 10 Mq black hole at unit Eddington rate M-Edd- The value of £ and A for which the calculation is performed 
are taken from the region PTR of Fig. 1. The scheme for the calculation is the following: 
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Fig. 3 Correlation among the QPO frequencies and the mass outflow rates for coupled accretion-wind system in <E>i(r). 

fnorm — ~ is 

the normalized dimcnsionlcss QPO frequency which is basically the ratio of the theoretically calculated QPO frequencies to the highest 
possible observed QPO frequency for any particular astrophysical source, see text for detail. 



First we calculate shock location r s h, shock compression ratio R CO m V and consequently the QPO frequency vqpo for the fixed 
lowest possible £ with the range of A for which the shock forms. We also calculate Rm for this set of parameters. We then 
repeat this calculation for higher values of energy till all £ and A covering the region PTR are exhausted. 
Fig. 3 is divided into four blocks. Each individual curve in any block represents the variation of fnorm on Rm resulted from 
the variation of A by keeping £ fixed, and different curves are drawn for different £ with A£ = 10 -3 . As stated earlier, 
£ and A covering the region PTR are used to draw Fig. 3. The whole range of specific energy of the shocked accretion 
material for which computation is done, was divided into four parts. The top left block represents results from very low energy 
accretion, top right for medium energy accretion, bottom left for moderately high energy accretion and bottom right is for 
very high energy accretion. One can observe that for very low energy accretion, normalized QPO frequency non-linearly and 
monotonically anti-correlates with the mass outflow rate and the calculated vqpo is also relatively low. As the specific energy 
of the accreting material increases, a peak appears in the correlation curve and the monotonic behavior is not maintained 
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Fig. 4: The outflow rate plotted against normalized QPO frequency for GRS 1915+105, as calculated from the observed radio flux and 
the measured QPO frequency. 



anymore. As £ becomes higher, the normalized QPO frequency first nonlinearly correlates with the mass outflow rate up to a 
certain value of vqpo (or with f„ or m) for which the mass loss rate becomes maximum, then the anti-correlation phase starts 
like low energy accretion cases. 

It is to be noted here that both fnorm as well as R m depend on initial boundary conditions and are quite sensitive to the 
shock location r B h- while f n0 rm non- linearly and monotonically anti-correlates with r 3 h, one can show that R m represents a 
distinctive non-monotonic variation with r s u- When shock location is small, R m nonlinearly correlates with r s h and attains 
a peak a particular value of r sh , say at r p ah . For r sh > r p ah (if available for a set of initial boundary condition), Rm starts 
anti-correlating with r s h- These informations might be used in the following way to infer whether our calculation is model 
dependent or not. 

For the same set of initial boundary conditions, the shock location distinctively differs for various $j used here. This is because 
of the fact that the region of parameter space (spanned by £ , A and 7) responsible for shock formation is quite different for 
almost all four $i (see Fig. 4 and Fig. 7 of D02). Hence for a particular value of [Pi], one obtains numerically different value of 
fnorm and Rm for all $i in general. This information, by any means, does not imply that our calculation is model dependent 
T because though the numerical values may differ, the overall fnorm — Rm (anti) correlation profile does remain the same for 
all $i. This is because of the fact that fnorm — r s h anti-correlation and R — m — r s h (anti) correlation profile remains unaltered 
for all $iS, only the numerical values of the related quantities change with the choice of potential. Thus we may, perhaps, 
conclude that our calculation is more or less 'generalized' in nature to include all existing pseudo-Schwarzschild potentials. 



t By the term 'model (in) dependent; we mean whether the basic physics remains unaltered or drastically changes with a specific choice 
of *j. 
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4 COMPARISON WITH OBSERVATIONS 

The Galactic micro-quasar GRS 1915+105 shows several properties which are indicative of the existence of shock oscillations in 
terms of QPOs and collimated jets in terms of radio emission (see Belloni 2002 for a general review on this particular source). 
Since it has showed intense and highly variable X-ray and radio emission characteristics, it has been extensively studied and 
it is an ideal source to test our concepts on shock oscillations and outflows. 

GRS 1915+105 shows diverse variability characteristics (see Belloni et al. 2000a for a classification) and the inter- 
relationship between these variability classes and the radio emission has been examined in detail (Naik & Rao 2000; Belloni 
et al. 2000b; Muno et al. 2001; Rao et al. 2000). Highly variable steep spectrum radio emission are invariably associated with 
superluminally moving ejecta and these are causally connected to a disturbed accretion disc as observed in X-rays as some 
sort of "inner disk oscillations" . These events could be initiated by some catastrophic events like the "magnetic rubber band" 
effect (Nandi et al. 2001). Steady flat spectrum radio emission, on the other hand, is always associated with the hard X-ray 
state of GRS 1915+105 and one of the identifying feature of such hard states is the existence of an ubiquitous 0.5 - 10 Hz 
QPO (Muno et al. 2001). It would be quite instructive to compare our results obtained in the previous sections to the outflow 
properties (as observed in radio) and shock oscillation properties (measured as the 0.5 - 10 Hz QPO) of GRS 1915+105. 

When the outflow is confined to the outer sonic point in the wind, the material can get collected and after some time 
can cause a catastrophic Compton cooling. This concept has been successfully used by CM to explain the repeated variation 
between high and low states in GRS 1915+105 and a relation between the QPO frequency and the 'low-state' (or 'off' state) 
duration has been derived and compared with the observations. Based on this concept, we can relate the predicted QPO 
frequency to the observed one by scaling to the maximum observed QPO frequency of 67 Hz (Morgan et al. 1997). 

To convert the observed radio luminosity to outflow rate, we adapt the method given in Fender & Pooley (2000). We 
assume that the radio spectrum has a flat spectrum from 1 GHz to 1.4x10 Hz (corresponding to mid-IR region of 2.3/im) 
and get a conversion factor of observed radio flux to jet luminosity as 1.35x 10 35 erg s _1 per mjy (for a distance of 12.5 kpc). 
The jet luminosity, Lj, can be converted to jet power, Pj by 

Pj^Ljn^F^i) (13) 

where r] is the radiative efficiency and F(T,i) is the correction factor for bulk motion (see Fender & Pooley 2000). rj~ 1 F(F,i) 
has a typical value of 10. Further assuming that the jet emitting region has a typical volume of 10 40 cm 3 and a magnetic field 
of 100 G, we can convert the jet power to outflow rate, getting a conversion factor of 1.4 e x 10~ 8 M yr _1 per 1 mjy for 
GRS 1915+105, where e is the fraction of particles accelerated to relativistic energies. In our present work, we have calculated 
the total outflow and further acceleration, presumably in some shock regions in the jet, would be required to accelerate these 
particles to relativistic energies. Since our purpose in the present work is only to get an estimate of the relative variation of 
outflow rate to shock oscillations, we have taken a conservative value of 0.1 for e. 

Using these conversion factors, we have plotted the observed outflow rate as a function of normalized QPO frequency, 
inorm-obs in Figure 4, based on a compilation by Muno et al. (2001). It should be noted that there is a general agreement with 
the trend of observations as predicted by the theory (see Figure 3) in the sense that the outflow rate shows an increase with 
tnorm at low frequencies, reaches a maximum, and then slowly declines. The agreement, however, is particularly close for the 
case of low specific energy (Figure 3, top left). The QPOs seen in GRS 1915+105 pertains to the steady low-hard states where 
the accreted matter can steadily loose its energy till it reaches the CENBOL and the specific energy can be very low for the 
accreting matter close to the black hole. The system can go into a high specific energy accretion mode at certain occasions 
(for example during the infrequent times when the high frequency - 67 Hz - QPOs are seen). It would be interesting to steady 
the behavior of QPO frequency and the outflow rate during such modes and we speculate that the variations would look 
more like the high specific energy behavior (Figure 3, bottom right). Also, there are possibilities that the outflowing matter 
(the ejecta) is essentially dominated by low binding energy (because they are supposed to be escaped towards infinity), hence 
such co-rrelation in stronger for low energy accretion. We must, however, caution that there are several uncertainties in the 
derivation of the outflow rate from the observed radio luminosity, but, the connection between the radio emission and the 
properties of the accretion disc seems to be on the line predicted by the Advective Flow formalism, described in detail in this 
paper. A detailed simultaneous observations and a calculation of the X-ray spectra using proper hydrodynamic formalisms 
would be required to further refine the model predictions. 

It is worth pointing out that the total amount of outflow observed in GRS 1915+105 during the 'baby-jet' episodes as well 
as during the super-luminal jet emission are also consistent with calculations presented here. Multi- wavelength observations of 
GRS 1915+105 (Mirabel et al. 1998) showed evidence for a blob of matter being ejected out, which are called the 'baby-jets', 
and it is estimated that > 10 19 g of matter is ejected during these events. It was pointed out by Belloni et al. (2000a) that 
these events are always preceded by hard dips characterized by ~3 Hz QPOs and they hypothesis that the matter for jet 
emission is ejected during these hard dips lasting for ~500 s. Hence, a mass outflow of ~10~ 10 Mq yr _1 must be participating 
in the relativistically accelerated particles, which is about 1% of the outflow estimated in our calculations. GRS 1915+105 
also shows superluminal jet emission with an estimated total mass of 10 23 g (Mirabel & Rodriguez 1999; Fender et al. 1999). 
These events could be either due to a series of dips preceded by low hard states with QPO frequency ~ 1 Hz (Naik et al. 
2001) or due to the accumulation of matter during an enhanced low state lasting for a few days (Dhawan et al. 2000). In 
either case, there is evidence that the ejected matter is accumulated during the low hard state with equivalent f„„™ of 0.01 
and lasting for ~50,000 s giving an effective outflow rate of 10 -8 Mq yr^ 1 , consistent with the calculation presented here. 
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5 CONCLUDING REMARKS 

In this paper we attempt to calculate the QPO frequency of the central accretor harbored by the galactic microquasar sources 
and to explore the dependence of the amount of baryonic load in the microquasar jets on these frequencies. Our calculation 
has been performed for polytropic shocked accretion onto galactic stellar mass black holes using a generalized post-Newtonian 
pseudo-Schwarzschild framework. However, one can show that shock formation is also possible when accretion is taken to be 
isothermal (Das, Pendharkar & Mitra 2003), and it is possible to calculate the normalized QPO frequencies for such accretion 
flow as well (Das 2003). In our future work, we would like to explore the correlation between vqpo and Rm for shocked 
isothermal inflow as mentioned above. 

Theory of accretion onto black holes is inherently a tough and intractable problem: first of all the nature of viscosity which is 
known to be turbulent even in mildly relativistic accretion flows is difficult to even formulate, let alone solve, in the relativistic 
accretion flows likely to be present in these system and the inner boundary condition is impossible to formulate and analytically 
handle in regions where general theory of relativity is known to predominate. The approach that we have taken in this work 
is the one in which general theory of relativity is handled in a set of pseudo-Schwarzschild potentials and the solutions are 
sought in regimes close to the black hole where the flow is shown to be supersonic and likely to be inviscid. Two fundamental 
problems are inherent in this approach namely: 

1) How are we sure that the pseudo-Schwarzschild potentials give correct description of general theory of relativity, 
and 

2) How do we extrapolate the results from inner regions of accretion disc where the accretion is likely to be inviscid to the 
outer regions where the flow is viscous and most of the observables like X-ray emission characteristics are measurable. 

We understand that such questions cannot be answered with mathematical uncertainty, but the present paper takes very 
important steps in that direction. By showing that shocks exist in all available pseudo-Schwarzschild potentials we have 
increased the confidence that shocks should indeed exist in the correct general relativistic formalism. By carefully selecting 
some important observables which are likely to be the manifestations of the conditions of inner accretion disk, we have 
attempted to make a semi-qualitative comparison with the theoretical predictions and have obtained very encouraging results. 
Though the fundamental problem of obtaining a complete solution for accretion onto black holes is not solved, very important 
insights are obtained in the present work which increases the confidence in the methodology and detailed prescriptions are 
given for comparison with observations. 
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